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Abstract 

A recipe is given for making materials with negative refraction in 
acoustics, i.e., materials in which the group velocity is directed oppo- 
site to the phase velocity. 

The recipe consists of injecting many small particles into a bounded 
domain, filled with a material whose refraction coefficient is known. 
The number of small particles to be injected per unit volume around 
any point x is calculated as well as the boundary impedances of the 
embedded particles. 
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1 Introduction 

Suppose that a bounded domain D C M 3 is filled with a material whose re- 
fraction coefficient n^(x) is known. Acoustic wave scattering problem consists 
of finding the solution to the equation 

[V 2 + k 2 n 2 (x)]u = in IR 3 , (1) 



u = e 



ika-x 



+ v: r(— ikv) = o(l), r — > oo. (2) 
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Here a G S is a given unit vector, the direction of the incident wave, S 2 is the 
unit sphere in M 3 , k > is the wave number, n^(x) = 1 in D' := M 3 \.D, and 
v is the scattered field. We assume that Imng(i) > 0. This corresponds to 
possible absorption in the material. The function u is the acoustic potential 
or pressure, and the time dependence factor e~ tujt is omitted. It is known 
that problem ([rj - (J2J) has a solution and this solution is unique. 

The problem we study in this paper is the following one: 
Can one create in D a material with negative refraction, that is, a material 
in which the group velocity is directed opposite to the phase velocity? 
If there is a wave-packet u = ^ k a{k)e l ^ k ' r ~' w ^ >t \ where \k — k\ < 5, \co(k) — 
uj(k)\ < 5, 5 > is a small number, a(k) is the complex amplitude, k is the 
wave vector around which the other wavevectors are grouped, and summation 
over k should be replaced by integration if vector k changes continuously, then 
the group velocity of this packet is defined by the formula: 

v g = V k u(k). 

Its phase velocity is defined by the formula: 

Vp ~ \k\ k ' 

where k° := Jh, \k\ is the length of the wave vector k, and Vfcu;(fc) is the 
gradient of the scalar function u(k). 

Suppose u{k) > 0. Then v p is directed along k. Vector v g = Vk^(k) is 
directed along —k if u = u;(|A;|) and w'(|A;|) < 0. Indeed, 

V fc ^(|A;|)=^(|A;|)A; , 

where the prime denotes the derivative with respect to the argument \k\. 

Therefore, one can produce the desired material with negative refraction 
if one can produce the material with 

u = u{\k\) and u'{\k\)<0. 

The notion of negative refraction in optics is discussed in detail in PQ, 
where the role of spatial dispersion in creating materials with negative re- 
fraction is emphasized. In [2] a theory of wave scattering by small bodies of 
arbitrary shapes is developed. On the basis of this theory it is proved in j3] 
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and in [3] that one can create material with a desired refraction coefficient 
n 2 (x) at a fixed frequency by embedding many small particles into the region 
D, filled with a material with known refraction coefficient Uq(x). 
Let us formulate one of the results from [I]. 

Suppose that the small particles are balls B m , 1 < m < M, M is the 
total number of small balls embedded in D, the radius of these balls is a, the 
distance d between two neighboring balls is of order 0(a 1//3 ), the boundary 
condition on the surface S m of the m— th ball is 

un = CmU on S m , 

where C, m := x m is the center of the m-th ball, and h(x) is an arbitrary 

given continuous function in D, with the property lmh(x) < 0. 
Assume that 



lim{aA^(A)}= [ N(x)dx 
J a 



for any subdomain A C D, where A/"(A) is the number of particles in A, and 
N(x) > is an arbitrary given continuous function in D. Let [b] denote the 
nearest integer to the numbe b > 0. 

Under these assumptions it is proved in J^j that embedding in every sub- 
cube A p C D the number [- J A N(x)dx] particles at the distances of order 
0(a 1//3 ) with boundary parameters ( m = h ^ Xm ^ ; creates new material in D 
with the refraction coefficient 



where 



n (x) = n (x) — k p(x), 

P(x) = ^p-N(x), xeD. (3) 
1 + h{x) 

Given p = p± + ip2, where p\ = Rep, P2 = Imp < 0, one finds (non-uniquely) 
hi(x) = Reh, h 2 (x) = lmh(x) < 0, and N(x) > 0, such that ([3]) holds. In 
[1] it was assumed that k was fixed, k = -, c = const is the wave speed in 
the exterior homogeneous region D' = R 3 \ D. 

The purpose of this paper is to point out that the method in J^j is valid in 
the case when p = p(x, u). This implies that this method allows one to create 
the material with the desired spatial dispersion. In particular, one can create 
a material with negative refraction in acoustics. 

Indeed, by the physical meaning of h one can prepare small particles with 
impedances ( = h( - x '^ ; depending on the frequency to in a desired way 
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Let us show that equation ([3]) can be solved non-uniquely with respect 
to the three functions: hi(x,u>), h 2 (x,oj) < 0, and N(x) > 0, given two 
arbitrary functions pi(x,u) and P2{x, uj) > 0. We emphasize that N(x) can 
be chosen so that it does not depend on uj. 

Lemma 1. For any pi(x,uj) and p2{x, uj) < 0, there exist three real-valued 
functions hi(x,u), h 2 (x,uj) < 0, and N(x) > 0, such that equation (j3J) is 
satisfied. The choice of these functions hi, h 2 < and N(x) > is non- 
unique. 

Proof. Our proof yields a method for finding a triple {hi, hi, N(x)} given the 
tuple {pi,p%}, and analytic formulas for the functions hi, h 2 , and N(x), see 
formulas (7), (13), (14) and (15). 
Denote 

R-= (pI+pI) 1/2 , 

and let ip be the argument of p = pi + ip 2 , so that 

Pi — Rcosip, p 2 = Rsinip; R — R(x, uj). (4) 
Equation (j3J) is equivalent to the following: 
47rJV(x)[/ii(l + /ii) + 



{l + hy + hi 

4tcN(x) h 2 



(l + hy + hi 

Let 



Rcosijj, (5) 
Rsxnip. (6) 



(1 + hi) 2 + h^ := r 2 , h 2 = rsimp, 1 + hi — r cos cp. (7) 
Equations (j3J) - (jSJ) can be written as: 

AttN(x)[1 coslp] = Rcosip, (8) 

47rN(x)^^ = R simp. (9) 



Let 



R \ p(x,u) if R(x, uj) > 0, 



AnN(x) 10 if R{x, uj) = 0. 



(10) 
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Then 

COS<p 

r 

sin ip 

r 

Given R and ip, we want to find r, ip and N(x) so that equations ( fTUl) - ( fT2 
hold. 

From dHJ) - (H2D one gets: 



p cos ^ — 1, (11) 
p siiaip. (12) 
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yp 2 — 2p cos -0 + 1 
If r is found by formula f[T3"j) . one finds 



(13) 



p sin ^ 1 — p cos if> . . 

sin ip = — - : , cos <p — — - . (14) 

\J p 2 — 2p cos if) -\- \ \J p 2 — 2p cos-^ + 1 

By definition (flQ]) p > 0. Choosing iV(a;) so that p < 1, one can make sure 
that 1 — 2p cosif) + p 2 > and cos<p > 0. The choice of N(x) > is made 
independently of u. 
One has: 

R(x, u) = 4:7tN(x) p(x, lu). (15) 

If p 2 < 0, which corresponds to some absorption in the material, then sinip < 
0, h 2 < 0, and sin <p < 0. If N(x) is sufficiently large, so that p < 1, then (fill) 
implies cos<p > and sin<p <0, so— |<(p<0. One can get the dispersion 
relation for u(k) from the equation 

u) 2 

— n 2 (x,u) = k 2 , (16) 
cr 

where c = const is the speed of the waves in the free space. 
From ffTBT) one gets 

u)n(x, lo) = \k\c. (17) 
Taking gradient with respect to k, one gets 

[n(x,u))+u)—]V k u) = ck°, k ° = J^- ( 18 ) 
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Therefore the group velocity 

Vg = V k UJ 

will be directed opposite to k° if 

u dn , , 

-^-<-l. 19 

In deriving ( Tl9l) we assume that n(x, u) is a real- valued function, or that 
the imaginary part of n(x,u) in D is negligible, so that when the wave is 
propagating through D its amplitude changes just a little. For this to happen 
it is sufficient to assume that L|Imn| <C 1, where L is the diameter of D. 

Since an arbitrary function n(x,u) can be created by the method of [I], 
one can ensure that inequality (fl9|) holds. This implies that the created 
material has negative refraction in the sense that the group velocity vector 
is directed opposite to the wave vector k. □ 

Example. If n(x,u) = (1 + ecu 2 )" 1 , where c > is a constant, then 
inequality (fT9l) holds if oJ 1 > 1. Suppose that the frequency interval is 
[uj m i n ,uj max \. Then inequality <HM holds for all frequencies from this interval 
if c> 

UT - 

771 % 71 

Let us give a summary of our method for creating material with the desired 
refraction coefficient: 

Step 1. Given the desired refraction coefficient n 2 (x,u>) in D, and 
the original refraction coefficient 71q(x,u), calculate the function p(x,u) = 
k 2 [nl(x,uj) — n 2 (x, u)]. 

Step 2. Given p(x,u), solve eq. (3) for h(x,u) = hi(x,ui) + ih 2 (x, u), 
h 2 (x,u) < 0, and N(x) > 0. 

There are many solutions: the three real- valued functions hi,hi < 0, and 
N(x) > are to be found from the two real-valued functions p± and p 2 < 0. 

In this step the solutions can be calculated analytically. 

Step 3. After h(x,cu) and N(x) are found in Step 2, one partitions the 
domain D into a union of small non- intersecting cubes Qj, centered at the 
points Xj, and embeds in each Qj the number Uj := [/ N(x)dx/a] of small 
balls of radius a, each having the boundary parameter Q = h(xj,uj) / a, where 
the symbol [b] denotes the closest integer to the number b > 0. 
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